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We apply the results of singularity analysis to the isotropic cosmological models in general relativity and string theory 
with a variety of matter terms. For some of these models the standard Painleve test is sufficient to demonstrate 
integrability or nonintegrability in the sense of Painleve. For others of these models it is necessary to use a less 
algorithmic procedure. 



1. Introduction 

In the modelling of various phenomena in relativistic 
cosmology the final result of the modelling process is a 
system of ordinary differential equations, whose solu- 
tion is by no means obvious. Generally these systems 
are of the first order and autonomous, so that, if the 
system is two-dimensional, one is assured of the exis- 
tence of solutions under mild conditions on the terms 
in the equations. For systems of dimension three or 
more the question of integrability or nonintegrability 
is extended by the possibility of chaotic behaviour in 
the general solution (for a recent review of the relevant 
dynamical systems methods in cosmology see p3)). 

One method to determine the integrability of a sys- 
tem is by a performance of the so-called singularity, 
or Painleve, analysis in an effort to examine whether 
or not there exists a Laurent expansion of the solution 
about a movable pole which contains the number of ar- 
bitrary constants necessary for a general solution. Any 
other singularities are not permitted except in the case 
of branch point singularities which give rise to what is 
called the weak Painleve property. A system which 
is integrable in the sense of Painleve has its general 
solution analytic except at the polelike singularity. 

The singularity analysis of ordinary differential 
equations goes back to the end of the nineteenth cen- 
tury to the pioneering works of Kowalevskaya [p^ , 
Painleve ||, ||, |2| and others |l|, § and has 
enjoyed a resurgence of interest in the last thirty years 
because of the growing attention paid to nonlinear 
equations, both partial and ordinary. The practical use 
of the singularity analysis has been greatly enhanced 
by the development of the ARS algorithm [Q, ||, ^ , the 
simplifications proposed by Kruskal ^ and the re- 
view by Ramani, Grammaticos and Bountis p6| . This 
practical approach has at times led to misinterpreta- 
tions at the hands of practitioners. The works of Conte 
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[H H have pointed the way to a rigorous application 
of the singularity analysis. 

In this paper, which is the fourth in a series ||l9, 20 
^ devoted to the investigation of connections between 
the three main topics in dynamics namely, symmetry, 
singularities and integrability, the topics which super- 
ficially are unrelated in their mathematics and yet are 
intimately intertwined, we apply the techniques of sin- 
gularity analysis to a number of model systems which 
have been derived for different problems in cosmol- 
ogy. The problems which we consider represent three 
classes, namely, those of fluid models in general rela- 
tivity, scalar field models coupled to fluid cosmologies 
again in general relativity and lastly isotropic string 
cosmologies. (For more information about the physical 
characteristics of the models treated below the reader 
is urged to consult the cited papers.) More specifically, 
we consider the following systems: 

1. The one-fluid FRW model in general relativity 
M (P- 122): 



H = -{l + q)H, 
h = -2q{l-n), 



g=i(37-2)f7 



(1.1) 



the latter of which is, of course, equivalent to the scalar 
equation 



h = (2-37)(l-f7)0. 



(1.2) 



Here, H is the Hubble variable, the density param- 
eter, q the deceleration parameter and derivatives are 
taken with respect to a dimensionless time parameter 
r defined so that dt/dr = 1/H, t being the usual 
time variable, giving the space sections of constant 
curvature in these models. As usual, the parameter 
7 € [0, 2] describes a barotropic fiuid. Although this is 
the simplest cosmological system, it is included here for 
completeness and comparison with more sophisticated 
ones introduced below. 

2. The two- fiuid FRW model in general relativity 
M (p. 126): 



ri = — i(fe — x) cos 2r2 cos 51, 
X = (1 -X^)smfl, 



(1.3) 
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where the so-called compactified density parameter 
fl G [— jTT, ^tt] and the transition variable x G 1] 
defined in the above reference are used to describe 
which fluid is dominant dynamically. Here b > —1. 

3. The flat one-fluid FRW space-time with a scalar 
field (f) with an exponential potential in general relativ- 
ity [pO| which reduces to the two-dimensional system 



T 2 



X = -3a; + 3(1- ^7) 



where differentiation is again with respect to a dimen- 
sionless time variable. Here A is a positive constant 
appearing in the scalar field potential and 7 S (0, 2) , 
that is, we exclude the "extreme" cases of stiff matter 
and scalar field coming from the fluid contribution. 

4. General-relativistic one-fluid FRW models with 
n scalar fields (f>i , i — 1, . . . ,n with exponential poten- 
tials described in dimensionless variables (fi, $, ^) 
by the system (this strictly excludes positively curved 
models wherein the variables are not defined) 



Z= 1,. 

i = 1,, 



where 



n n 

g=i(37-2)f^ + 25]*,?-^$,^ 



(1.5) 



(1.6) 



i=l i=l 



The variables (<&, 5*) are defined in terms of the po- 
tential and kinetic energies of the scalar fields and 
7G[0,2]. 

5. The four-dimensional flat string FRW model 
with negative central charge deficit described in the 
compactifying variables (^, r]) by the system [B2| 



(1.7) 



6. The ten-dimensional flat string FRW model in 
the RR sector and with a positive cosmological con- 
stant described by the system 



X 

V 
z 



{x + V3)(l - - y - z) + \z{x - \/3), 
2y{{x ^ %/3)(l - - y - z) + \z{x - %/3)], 
2z\[x ^\fZ){\- x^ -y- z)^ \z{x - Vs)], 

(1- 



for which the dimensionally reduced (four-dimensional) 
model is obtained by placing y = in (1.8). 

Another approach to determining the integrability 
of sets of differential equations is the use of Noether's 
theorem for Lagrangian systems and the Lie theory of 



extended groups for differential equations in general. 
For the systems which we study here the Lie symmetry 
approach is not generally successful because we need to 
find generalised or nonlocal symmetries to supplement 
the point symmetries which are the easiest to obtain. 

For the benefit of the reader we describe in Sec. 2 
the main points of the Painleve process (for a more 
detailed albeit elementary introduction we refer the 
reader to ||2^). This is then used in Sec. 3 to study 
the singularity features of our cosmological systems. 
Conclusions are given in Sec. 4. 

2. Methodology of Painleve analysis 

The singularity analysis which lies at the basis of the 
Painleve test as systematised in the ARS algorithm [1- 
3] is a specific form of a more general analysis which ex- 
amines a set of differential equations for leading order 
behaviour and next to leading order behaviour |l^, |2j . 
The latter analysis does not take into consideration the 
necessity for a Laurent expansion about a polelike sin- 
gularity (or rational branch point) which is required 
for the Painleve test. The essence of the Painleve test 
(for a set of ordinary differential equations which is the 
only type which we consider here) is that the solution 
of an n-dimensional set of equations 



x-i — fi (i, x) , 



(2.1) 



where the functions /; are rational in the dependent 
variables and algebraic in the independent variable, 
can be written as either 



Xi{t) = > UiT 



00 



in the case of a Right Painleve Series or 



00 

J=0 



-Pi-gj 



(2.2) 



(2.3) 



in the case of a Left Painleve Series |T^, where 
T — t — to and to is arbitrary, the exponents pi are 
strictly positive integers (rational numbers in the case 
of the so-called weak Painleve test), the parameter q is 
a (positive) integer (respectively rational number) and 
in the coefficients there are n — 1 arbitrary constants 
which, together with to , give the required number of 
arbitrary constants for the general solution of the sys- 
tem (^). 

The ARS algorithm provides a mechanistic proce- 
dure for the determination of the leading order be- 
haviour and the resonances which are where the ar- 
bitrary constants of integration arise. The first step is 
to determine the leading order behaviour by means of 
the substitution 

X, = a,T-P^ (2.4) 
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into the system (|2.lD and to assemble all possible pat- 
terns of values of the exponents pi,i = l,n from the 
dominant terms of the system (2.1). For each of these 
possible patterns the coefficients at are determined 
and the substitution 



(2.5) 



is made to determine the "resonances" r at which 
the arbitrary constants, /^i , are introduced. The res- 
onances and the arbitrary constants are determined 
from the leading order behaviour of the system ( p.l| ) 
by a linearisation process. The final step for each par- 
ticular pattern of leading order behaviour is to substi- 
tute a series, truncated at the highest resonance, to 
ensure that there is compatability at the resonances. 
If, for a particular pattern of leading order behaviour, 
these conditions are satisfied and the series contains 
n — \ arbitrary constants (the nth comes from t^), 
that particular pattern is said to pass the Painleve test. 
If this is true for all possible patterns of leading order 
behaviour, the set of equations is said to possess the 
Painleve property and to be integrable. We note that 
the concept of Painleve integrability means the posses- 
sion of a Laurent series, i.e. , the solution is analytic 
except at the polelike (branch point in the case of the 
weak property) movable singularities. 

The algorithm described in the previous paragraph 
is not always possible or easy to implement. The first 
problem arises when not all of the exponents obtained 
in the leading order analysis are strictly positive. In the 
case that some or all of the exponents are strictly neg- 
ative one can perform a homeographic transformation, 
which preserves the Painleve property, on the affected 
variables so that the exponents will now be strictly pos- 
itive. The test can then be continued algorithmically. 
There is no such consolation in the case of zero expo- 
nents in the leading order behaviour. The algorithm 
fails. In this case it is necessary to make a series sub- 
stitution commencing at the leading order behaviour 
to establish whether or not a sufficient number of arbi- 
trary constants is introduced. If these constants enter 
at early terms in the series, the task is not too difficult 
with one of the symbolic manipulation codes. If this 
is not the case, the task can become impossible for a 
system of moderate complexity. 

Another possibility which can occur and create a 
problem of interpretation is that the number of arbi- 
trary constants introduced at the resonances is insuf- 
ficient to produce the general solution for the system. 
This possibility was already observed by Ince, citing 
from a paper by Chazy written in 1909, in 1927 
(p. 355) and treated in some more detail in pO| . Ince 
described such a solution as a "singular" solution. This 
is in accord with the usual usage of the word in the case 
of first order equations. However, it is a little unfortu- 
nate that the word "singular" has to perform two dis- 
tinct functions in the description of the integrability of 



one system. In an attempt to separate the two uses of 
the word, Cotsakis and Leach jl^ described this type 
of solution as a "partial" solution confined to a sub- 
manifold of the space of initial conditions, but this has 
not met universal acceptance. One is reluctant to use 
"particular" solution because of the commonness of the 
term in the solution of nonhomogeneous linear equa- 
tions. In Greek there are three words, " a^^w/iaAjy" , 
"eiJiKry" and " tJidCoucra" meaning "singular", "par- 
ticular" and "peculiar" respectively, which in Greek 
are used to describe these different types of solution. 

In this paper we propose to use the adjective "pe- 
culiar" to describe a solution of the type given by Ince. 
As far as we are concerned here, the question is whether 
or not the existence of such solutions violates the re- 
quirement that all possible solutions pass the Painleve 
test, which, one must recall, requires the existence of 
the correct number of arbitrary constants in the Lau- 
rent expansion for each pattern of leading order be- 
haviour. According to Tabor |2^ (p. 330), this is 
necessary for the possession of the Painleve property. 
However, a recent paper |l^ has demonstrated inte- 
grability in the case that of the two possible patterns 
of leading order behaviour one satisfied the Painleve 
test and the other was lacking one arbitrary constant. 
This demonstration was for a single example and one 
would want a sounder basis for making a definite claim 
about integrability in cases for which the requirements 
of the Painleve property were only partially satisfied 
for some patterns of leading order behaviour. 

Finally we recall that the possession of the Painleve 
Property is representation dependent and is definitely 
preserved only under a homeographic transformation. 
However, there are times when we can introduce a 
transformation with profit such as in the system de- 
scribed by (1.3). The nature of the integrability of the 
original system in comparison with that of the trans- 
formed system will be tempered by the nature of the 
particular transformation. 

3. Applications of Painleve analysis to 
models 



3.1. One-fluid FRW model, ( |l.l| ) 



In the case of (1.1) it is useful to introduce the rescaling 
transformation 



(3.1) 



H — >x, n — 

to give the system 
X = -{(T + y)x, 

y = -2y(l-y), (3.2) 

where a = ^ (37 — 2) . When we make the usual sub- 
stitution for the leading order behaviour, we obtain 

Pqr"-^ = 2/3V29. (3.3) 
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From the second Eq. (B.2) it is evident that q = — 1 
and /9 = — i . From the exponents of the first Eq. (3.2) 



the value of p is arbitrary. However, equahty of the 
coefficients of the leading order powers imposes the re- 
quirement that p — ^ . Although the value of p is not 
strictly negative, we can proceed with the Painleve test 
without a further transformation since it is a noninte- 
gral rational number. We note that the value of a is 
unspecified. 

To obtain the resonances we make the substitution 



X 

y 



ar^/^ + pT-+^/\ 



r-l 



to obtain the condition 



a 







(3.4) 



(3.5) 



that there be a nontrivial solution. The condition (3.5) 
gives the resonances r = —1,0 and so the system (3.2) 
passes the Painleve test. Hence the system (1.1) pos- 
sesses the Painleve property. The first few terms of the 
Laurent expansion are 



x{t) 



.3/2 



Ma- 



5/2 



6 



(3.6) 



3.2. Tv^fo-fluid model, 



The system (^j) is not in a suitable form for applying 
the Painleve analysis. We have two choices for the in- 
troduction of new variables. Interestingly both changes 
of variables produce essentially the same results. In the 
first instance we introduce the new variables 



X 



and 



y — cosecfi, 



so that the system (O) becomes 



xy 



= \{b^x){y^^l){y'-2). 



(3.7) 



(3.8) 



The usual leading term analysis gives the exponents 
p = — 1 and q = 0, so that the ARS logarithm is not 
applicable. 

We do not continue with the analysis of the system 
in (3.8) since, as we noted above, the second change of 



variables produces essentially the same results. 

The second transformation, which we consider in 
detail here, has the advantage of being one-to-one and 
continuous over the defined intervals of the original 
variables. We set 



X — X a-iid y = sin 
to obtain the system 
i = (1 - x^)y 

y = -^(b-x){l-2y')il-y'). 



(3.9) 



(3.10) 



The same exponents are obtained as for (3.8). We 
make the Ansatz 



i=0 



(3.11) 



and substitute this into the system (3.10) to obtain the 
pair of relations 



(3.12) 



ib,,T 



i-l 



-i&(l - 3boibjT'+^ + 2hbjbkbiT'+^+^+^) 
+ a,(r'-i-36,6fer*+^+^-i 
+ 26,6fc&i6,^T'+^+'=+'+™-i), (3.13) 

from which we are able to deduce the first few terms of 
the expansions for x and y. From the first two terms 
of ( 3.12 ) we obtain 



apfeo — 1 and oi = —^a^bi 
and from the first term of ( 3.13| ) 
ao (1 - ihl + 2hl) = 0. 



(3.14) 



(3.15) 



Since 7^ 0, it follows from (3.15) that — 1 



In either case the second term of ( 3.13| ) gives bi — 
and consequently ai — Q. The third term of ( 3.13| ) 
reduces to 62 = ^§^2 , so that, for the first possibility 
for the value of &o i ^2 is arbitrary and, for the second 
possibility, zero. In fact for the second possibility all 
subsequent coefficients hi are zero and we obtain the 
solution 



x{t) 



1 

bo, 



bpT 
3 



(bor) 
45 



(3.16) 



where 6q = 1/2, which is certainly a peculiar solution. 

In the case that 6g = 1 we obtain a more standard 
solution. We have that 62 is arbitrary and this pro- 
vides us with the second arbitrary constant required 
for a general solution of the original system. We have 



aoT-^ + i (foo - 262) T + lbbob2T^ 
bo + b2T^ + . . . , 



(3.17) 



where the coefficients ao and bo have been given above. 

In the expansions we have obtained, only the sec- 
ond one has the required number of arbitrary con- 
stants, and we cannot conclude that the system ( |l.3| ) 
is integrable in the sense of Painleve. However, we 
do note that it is possible to obtain a first integral of 
the original system (^]^) and reduce the solution to a 
rather complicated quadrature. In fact the system can 
be written in terms of a Lagrangian and is Hamilto- 
nian, so that the existence of the first integral immedi- 
ately guarantees integrability in the sense of Liouville. 



Painleve Integrability of Isotropic Cosmologies 



5 



3.3. Flat FRW with one fluid and an 



exponential potential, (1.4) 



Before we begin the singularity analysis of (1.4), it is 
appropriate to simplify the system by the rescaling 

AV3 



V2(2-7) ' 
Py with = 



3A2 



27(2-7)' 



(3.18) 



in which the sign of (3 may be taken as positive without 
loss of generality. The system ( |l.4| ) now has the simpler 
appearance 



xy 



y = ~xy + x^y - y^, 



(3.19) 



where ^ = 8 (l - ^7) / (3A2) and B = AV (7(2 - 7)) . 

We determine the leading order behaviour of the 
system ( |3.19| ) to be x{t) ^ cxt-'^I'^ and y(r) = /Jr-^/^ 
with the constraint c? — (? — —1/2- We find that the 
resonances are at r = — 1, 0, where the first resonance 
is generic and the second indicates that one of the co- 
efficients of the leading order behaviour is arbitrary, 
which is in accordance with the above constraint. If 
we make the substitutions 



.(»-l)/2 



x{t) = J2' 

1=0 



(3.20) 



we find that the first few terms of the expansion are 
given by 



ao 
ai 

a2 



"0 
61 
b2 



arbitrary 

B + l{2 + 5B)al 

iao(3S(l + 5S)- 

+28(1 + S)(2 + 

+80(1 



+ 1(1 
1-8(14 
5B)a^ 
- 6A{1 



VB)4 
25(4 -f 



bB))al 



|ao(-l 



+ 2B + A{\ + B)al) 
l + B)B 
+2(-2 -3A + 2B + 16B^)al 

+ B){2 + 23B)a^ 
+80(1 + Bfal). 



(3.21) 



We conclude that the system (1.4) is integrable in the 
sense of Painleve. 



3.4. Exp onential potential with one fluid 

We consider the case in which n = 1. When we make 
the substitution for q in the system ( |l.5| ), we obtain 

ii = 1(37 - 2)*r2 + 2*3 - - 2* - \/iA"$^ 
$ = i(37-2)$r2 + 2*2$_$3^^_ 

17 = (37 - 2)n^ + - 2<P^n - (37 - 2)n. 

(3.22) 

We find that the leading order behaviour is given by 
«'(t) = ar-^/^, $(r) = br'^^^, 17(t) = cr tJ,23) 

subject to the constraint that 

(37- 2)c + 4a2 - 26^ = -1. (3.24) 

The analysis of the dominant terms for the resonances 
is facilitated by the substitutions 



4*^ ^; = -2$^ 



{3j-2)n. 



(3.25) 



(Note that this transformation does not preserve the 
Painleve property. However, it is satisfactory for the 
purposes of this immediate analysis.) 



The dominant terms in the system (3.22) are now 

X — x{x + y + z) 
y = y{x + y + z) 

z = z{x + y + z). (3.26) 

The leading order behaviour of ( 3.26| ) is given by 
.-1 



X = ar 



7t"\ 



(3.27) 



subject to the constraint a+/3+7 = — 1 . (The constant 
7 in (B.27) is not to be confused with the physical 
constant in the original system.) We determine the 
resonances by substituting into (B.26) 



X = OLT 

y 

z 



-/IT , 

(3t^^+i^t'~\ 



7t"^ + pr*~^ 



to obtain the linearised system 



(3.28) 



0(3.29) 



which has a nontrivial solution if s = —1,0(2). Thus 
we see that there is a double zero resonance which is 
consistent with the constraint. These results pass over 
to the original system and, since two arbitrary con- 
stants enter at the leading order terms, the system 



s — a 


—a 


—a 






-a 


s-P 






V 


-7 


-7 


s — 7 




P 
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{ \i-22\ ) passes the Painleve test for this pattern of lead- 
ing order behaviour. 

We present the first few terms of the Laurent ex- 
pansion of the original system (3.26). 
1 

+ ai 



1 

a2T2 4 

1 

+ 6i +62T2 + 



* = aoT 2 

where 

ai = -\/6fc (1 + ial 



C2 



(3.30) 



02 



hi 



^ao[-8 + A-2{2 + A + 12k^) hi 



+al [(-4 + A + l^k%l) (-1 + 106^)] 
Vefcao&o (1 - 462) 



-2{2 + A + &k^{l + l8al))bl 
-|-36fc2 (1 + 20ao) 60] 



Co 



1 



-l-4a2-F2&2) 



Cl 



C2 



-^[{l^Aal-2bl){-A 

-IQal -I- 4Aag - 2(2 + A + mk^aDhl 
-^12fc2(3-f 76a^)6;5)], (3.31) 

in which we have taken ao and &o as two arbitrary 
constants and expressed the third, cg , in terms of them 
using the constraint ( ^.24 ) . 



We note that the system ( 3.26 ) is easily integrated 
to give the closed-form solution 

1 



x{t) 

y{t) 

z{t) 



(1+Ci+C2)(i-C3) 

C2 

(1+Cl+C2)(i-C3) 
Cl 



(3.32) 



(1 + ci -I-C2) (i- C3) 

Hence the leading order analysis yields the exact so- 
lution of the system consisting of the dominant terms 
only. 

Apart from the leading order behaviour given in 
( ^.27 ), there are other possibilities. Denoting the lead- 
ing order powers of x , y and z hy p, q and r , we have 
the following possibilities apart from that in ( ^.27 ) 
listed in the table: 



p 


9 


r 


-1 


-1 


> 


> 


-1 


-1 


-1 


> 


-1 


-1 


> 


> 


> 


-1 


> 


> 


> 


-1 



In each of these cases it is necessary to make a series 
substitution since the ARS algorithm for the applica- 
tion of the Painleve test is no longer appropriate. We 
summarise the results: 

1. We substitute the series 



i=Q 



i=0 



(3.33) 



to find that ao is arbitrary and cq — ci — C2 — ■ ■ ■ — , 
so that we have a peculiar solution with two arbitrary 
constants (ao and to) and one of the functions, c(t) , 
having only the trivial solution. 

2. The series substituted are now 



1=0 



1=0 



-(■'-2)/2 



(3.34) 



We find that uq and cq are arbitrary and that bo = 
^(1 + Aco)/2, so that the series solutions do contain 
three arbitrary constants when to is included. 
3. In this case we put 



a/2 



i=0 



j=0 



-(■'-2)/2 



(3.35) 



1=0 



We distinguish three subcases. 

Subcase (a): the coefficients ao = ai = a2 = . . . = 
and 60 = ^1 = ^2 = • • • = and cq — ^1/A indi- 
cate that we have trivial solutions for ^ and $ and a 
nontrivial series with arbitrary constant to for Q . 

Subcase (b): we find that ao is arbitrary and bo = 
bi ~ b2 ~ ■ ■ ■ — , so that we have the trivial solution 
for $ and a two parameter solution ( ao and to ) for ^ 
and Q . 

Subcase (c): In this subcase ag = ^i, bo — bi = 
62 = • • ■ = and Co = Cl = C2 = . . . = , so that we 
have trivial solutions for (f> and Q and a one-parameter 
solution for ^E" . 

All three subcases present us with peculiar solu- 
tions in the form of series. 

4. With the substitution 



.1/2 



1=0 



i=0 



i7 = E< 



i/2 



(3.36) 



we find that ao = 0, so that there is no singular be- 
haviour in any of the series representations of the three 
functions. However, the constants ai , bo and co, the 
leading terms of each of the series, are arbitrary, and 
so we obtain a three-parameter representation of the 
solution. 
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5. When we make the substitution 

i=0 1=0 

n = J2c,T^/^ (3.37) 



i=0 

we find two subcases. 

Subcase (a): the coefiicient = 0, which removes 
any possible singular behaviour, and the coeflicients 
ao , bi and cq are arbitrary, thereby providing a three- 
parameter series representation for , <i> and il . 

Subcase (b): the coefficient bo ~ 2^^/^. All coeffi- 
cients Ci are zero and all coefficients Oi and &i, z > 
are expressed in terms of the parameter A of the sys- 
tem. Thus we have a trivial solution for f2 and a one- 
parameter solution (the location of the singularity to 
of $ ) for the two functions ^' and $ . 

The first subcase presents a solution without a sin- 
gularity and the second one a peculiar solution with 
singularity. 

6. For the last case we substitute 



i=0 



i/2 



i=0 i=0 
.(»-2)/2 



(3.38) 



i=0 



We find that cq — ci — 0, thereby removing any pos- 
sible singular behaviour, and that ao , bo and C2 are 
arbitrary. Thus we have a three-parameter series rep- 
resentation of the solutions for 4" , $ and fl . 

Altogeth er t he results support the proposition that 
the system (L5) is integrable. 



3.5. Four-dimensional flat string FRW, (1.7) 



In the analysis of the leading order behaviour of the 
system (1.7) we find the leading order behaviour 



V 



(XT 



1. 



(3.39) 



In accordance with the constraint in (3.39) we find that 
the resonances are at r — —1,0. Since the second 
arbitrary constant enters at the leading order terms, 
the full system, (1.7), passes the "weak" Painleve test 
and so is integrable in the sense of Painleve. The first 
few terms of the Laurent expansion are 



1 



-l + 3\/3bg-H86 ^ ^ 
1564 



/3 



' 352rV3 
+ . . . 

_ 3(-l + 3V3b-] + 3b[S) 
^ t1/3 155;! ^ 



(3.40) 



(3.41) 



3.6. 10-dimensional flat string FRW, (|TT8| ) 

For the system (O) 



= (x-f V3)(l-a;^ 



22/[(l 



z) -f \z 



z ^ 2z 



(1 



z) + \z{x-^/?,), (3.42) 
(3.43) 

(3.44) 



- y - z) - \{\ - z - Vi:, 



we determine the leading order behaviour of this sys- 
tem to be x{t) — aT"^!"^ , y (r) — /3r^-^ and z(t) = 
7T^^ with the constraint 2c? + 2/3 + 7 = 1. We deter- 
mine the resonances by the substitution of 

y(T) = /3T-i+j.r^-i, 
Ar) - ' ' 
into (F 



7r + pr 
to obtain the linearised system 



(3.45) 



s + 2a2 


a q;/2 \ 






4a/3 


s + 2/3 /? 






4q7 


27 s + ^ I 







(3.46) 



which has a nontrivial solution if s = — 1, (2) . There 
is a double zero resonance which is consistent with the 
constraint and the system passes the Painleve test for 
this pattern of leading order behaviour. The first few 
terms of the Laurent expansion 

X (r) = aoT"i/2 j^ai+ oar^/^ + . . . 
y (r) = boT-^ + + 5^ + ^3^1/2 

rl/2 



z (r) = Cor ^ + ciT + 
are given by 

ai = il/Vt] [3 (-3 + 460) 



C2 + C3T 



(3.47) 



02 

Co 
Cl 

C2 



-2al{-ll + &al + &bo)] 
(ao/2) [-93 + ISQal + 46o (68 - 456o) 

+6a^(-77 + 60&o) 

+2al (188 + ibo (-107 + 306o))] 
- (4/V3) aobo (-5 + Qal + 660) 

(feo/3) [-63 + 168feo + 2(342a|^ - 546o 

+3a;5 (-199 + 22860) 

+al (290 + 36o (-217 + 1146o)))] 
1 - 2ao - 26o 

(2/V3) ao (-1 + 2a2 + 26o) (-13 

+ I2al + 1260) 
(-1/3) (-l + 2ag + 26o) [QMal 
+6ao (-229 + 2286o) 
+2ag (392 + 576o (-13 + 660)) 
-3(31 + 46o(-17 + 96o))]. 



(3.48) 



Note that the constraint is apparent in the first of 
the relations for the coefficients c, . 
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4. Discussion 

The examples discussed in this paper may prove useful 
when one is interested in having a theory to decide the 
question of what is the general singularity pattern of 
isotropic cosmologies in different gravity theories. This 
question, apart from its purely mathematical interest, 
is believed to be related to recent observations of the 
possible oscillatory nature of the universe, exemplified 
in an oscillatory behaviour in the Hubble parameter 

i- 

Among the models discussed in this paper, the 
two-fluid FRW model in general relativity presents 
the most interesting dynamical behaviour. All other 
models are either strictly integrable in the sense of 
possessing the strong Painleve property, or they have 
branch-point singularities indicating weak integrability 
in the sense of Painleve. 

The two-fluid model possesses the interesting prop- 
erty of what can be called singular envelopes, first dis- 
cussed by Ince and rediscovered in a different context 
in ijl^. That is, although the general solution is un- 
known, any possible nonintegrable or even chaotic be- 
haviour may be confined to a region of phase space 
enveloped by the peculiar solutions in the sense given 
at the end of the previous section. It is interesting to 
ask whether this property is rigid, that is, is it main- 
tained when two-fluid models are considered either in 
other gravity theories, or in more general Bianchi mod- 
els in general relativity. The former question is cur- 
rently under investigation and can be reformulated as 
a two-fluid plus a scalar field model in the Einstein 
frame. 

Another question is whether integrability in the 
case of scalar field models is maintained when one con- 
siders more general potentials. This is known not to 
be the case even for a general FRW model with a 
scalar field that has a simple quadratic potential 0. 
One would like to be able to relate the integrability 
properties of different cosmological spacetimes in the 
context of different gravity theories and matter fields 
in an effort to understand the significance of excep- 
tional non-integrable cases as opposed to generic, in- 
tegrable ones in the simple frame of isotropic models 
before one moves on to more difficult homogeneous but 
anisotropic case. Problems in this direction are cur- 
rently under investigation . 

Acknowledgement 

POLL thanks the Director of GEODYSYC, Dr S Cot- 
sakis, and the Department of Mathematics, University 
of the Aegean, for the use of their facilities during the 
time when this work was performed and acknowledges 
the continuing support of the National Research Foun- 
dation of South Africa and the University of Natal. 



References 

[1] M.J. Ablowitz, A. Ramani & H. Segur, Nonlinear evo- 
lution equations and ordinary differential equations of 
Painleve type Lett Nuovo Cimento 23, 333-338 (1978) 

[2] M.J. Ablowitz, A. Ramani & H. Segur, A connection 
between nonlinear evolution equations and ordinary 
differential equations of P-type. I J Math Phys 21, 
715-721 (1980) 

[3] M.J. Ablowitz, A. Ramani & H. Segur, A connection 
between nonlinear evolution equations and ordinary 
differential equations of P-type. II J Math Phys 21, 
1006-1015 (1980) 

[4] T.J. Broadhurst, R.S. Ellis, D.C. Koo and A.S. Szalay, 
Nature 343, 726-728 (1990) 

[5] J. Chazy, Sur les equations differentielles du troisieme 
ordre et d'ordre superieur dont I'integrale generale a 
ses points critiques fixes Acta Math 34, 317-385 (1911) 

[6] P.A. Clarkson & M.D. Kruskal, New similarity reduc- 
tions of the Boussinesq equation Nonlinearity 7, 975- 
1000 (1989) 

[7] A. A. Coley & R.J. van der Hoogen, Oscillatory be- 
haviour in cosmological models / J Theor Phys 36, 
2293 (1995) 

[8] R. Conte, Universal invariance properties of Painleve 
analysis and Backlund transformation in nonlinear 
partial differential equations Phys Lett A 134, 100-104 
(1988) 

[9] R. Conte, Singularities of differential equations and 
integrability, in Introduction to Methods of Com- 
plex Analysis and Geometry for Classical Mechanics 
and Nonlinear Waves Benest D and Froeschle C edd 
(Editions Frontieres, Gif-sur-Yvette) 49-143 (1994) 

[10] S. Cotsakis & P.G.L. Leach, Painleve analysis of the 
Mixmaster universe J Phys A: Math Gen 27, 1625- 
1631 (1994) 

[11] SW. Cotsakis, P.G.L. Leach & J. Miritzis, Painleve 
integrability of scalar field cosmologies (preprint: 
GEODYSYC, Department of Mathematics, University 
of the Aegean (2000)) 

[12] M.R. Feix, C. Geronimi, L. Cairo, P.G.L. Leach, R.L. 
Lemmer & S. Bouquet, On the singularity analysis 
of ordinary differential equations invariant under time 
translation and rescaling J Phys A: Math Gen 30, 
7437-7461 (1997) 

[13] B. Gambier, Sur les equations differentielles du second 
ordre et du premier degre dont I'integrale generale est 
a ses points critiques fixes Acta Math 33, 1-55 (1889) 

[14] R. Garnier, Sur les equations differentielles du 
troisieme ordre dont lintegrale generale est uniforme et 
sur une classe d'equations nouvelles d'ordre superieur 
dont I'integrale generale a ses points critiques fixes 
Ann Set Ecole Norm Sup XXIX, 1-126 (1912) 

[15] E.L. Ince, Ordinary Differential Equations (Long- 
mans, Green & Co, London (1927)) 

[16] M. Jimbo, M.D. Kruskal & T. Miwa, Painleve test 
for the self-dual Yang- Mills equation Phys Lett A 92, 
59-60 (1982) 



Painleve Integrability of Isotropic Cosmologies 



9 



[17] S. Kowalevski, Sur le probleme de la rotation d'un 
corps solide autour d'un point fixe Acta Math 12, 177- 
232 (1889) 

[18] P.G.L. Leach, Hierarchies of similarity symmetries and 
singularity analysis Dynamical Systems, Plasmas and 
Gravitation P G L Leach, S E Bouquet, J-L Rouet & 
E Fijkalow edd (Springer- Verlag, Heidelberg) 304-312 
(1999) 

[19] P.G.L. Leach, S. Cotsakis & G.P. Flessas, Symmetry, 
singularities and integrability in complex dynamics 1: 
The reduction problem J Nonlm Math Phys 7, 1-34 
(2000) 

[20] P.G.L. Leach, S. Cotsakis & G.P. Flessas, Symmetry, 
singularities and integrabily in complex dynamics II: 
Rescaling and time-translation in two-dimensional sys- 
tems J Math Anal Appl (to appear) 

[21] P.G.L. Leach, S. Moyo, S. Cotsakis & R.L. Lemmer, 
Symmetry, singularities and integrability in complex 
dynamics HI: Approximate symmetries and invariants 
J Math Anal Appl (to appear) 

[22] R.L. Lemmer & P..L. Leach, The Painleve test, hidden 
symmetries and the equation y" +yy' + ky^ = J Phys 
A: Math Gen 26, 5017-5024 (1993) 

[23] P. Painleve, Memoire sur les equations differentielles 
dont I'integrale generale est uniforme Bull Math Soc 
France 28, 201-261 (1900) 

[24] P. Painleve, Sur les equations differentielles du second 
ordre et d'ordre superieur dont I'integrale generale est 
uniforme Acta Math 25, 1-85 (1982) 

[25] P. Painleve, Sur les singularaties essentielles des 
equations differentielles d'ordre superieur CRAS 116, 
362-365 (1893) 

[26] A. Ramani, B. Grammaticos & T. Bountis, The 
Painleve property and singularity analysis of inte- 
grable and nonintegrable systems Phys Rep 180, 159- 
245 (1989) 

[27] M. Tabor, Ghaos and Integrability in Nonlinear Dy- 
namics (John Wiley & Sons, New York (1989)) 

[28] J. Wainwright, in General relativity, Proceedings of 
the Scottish Universities Summer School, 46 Hall J S 
& Pulham J R edd (lOP Publishing, Bristol (1995)) 

[29] A. A. Col ey, Dynamic al systems in cosmology 
(preprint: |gr-qc/991007^ 

[30] E.J. Copeland, A.R. Liddle & D. Wands, Exponential 
potentials and cosmological scaling solutions Phys Rev 
D 57, 4686 (1998) 

[31] A. A. Coley & R.J. van der Hoogen, The dynamics 
of multi-scalar fiel d cosmological m odels and assisted 
inflation (preprint: gr-qc/9911075 ) 

[32] A. P. Billyard, A.A. Coley & J.E. Lidsey, Quali- 
tativc analy sis of string cosmologies (preprint: ^r^ 
qc/9903095^ 

[33] J. Wainwright and G.F.R. Ellis, Dynamical Systems in 
Gosmology (Cambridge University Press, Cambridge, 
1997) 



